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ABSTRACT 

-A 

Analytical  and  experimental  data  is  presented  validating 
a  finite  segment  cable  model.  The  model  consists  of  a  series 
of  pin-connected  rigid  rods  which  may  have  different  lengths, 
diameters,  and  masses.  The  model  is  capable  of  simulating  large, 
three-dimensional  motion  of  flexible  cables.  Its  principal  areas 
of  application  are  expected  to  be  with  the  simulation  of  long, 
heavy,  towing  and  hoisting  cables. 

x  - 


INTRODUCTION 


Approximately  five  years  ago  we  introduced  a  three-dimensional 
dynamic  finite  segment  model  for  cables  and  chains  [1] .  *  The  model  consists 
of  a  series  of  pin-connected  rods  or  links.  The  model  is  developed  so  that 
arbitrary  external  forces  may  be  exerted  on  the  links.  The  dynamics  of 
the  system  is  then  determined  through  numerical  integration  of  the  governing 
equations  of  motion.  In  this  report  we  present  data  which  validates  ana¬ 
lytical  predictions  of  the  model. 

Figure  1.  shows  a  schematic  representation  of  the  model.  The  number 
of  links  N  is  arbitrary.  The  length,  diameter,  mass,  and  inertias  of  the 
individual  links  may  also  be  arbitrarily  chosen.  This  model  is  expected 
to  be  effective  in  studying  the  nonlinear,  three-dimensional,  dynamic  behavior 
of  long,  heavy  cables  [1,2] .  Particular  application  with  submerged  towing 
cables  has  also  been  suggested  [3] . 

The  attractive  features  of  the  model  are:  1)  the  arbitrary  dimensions 
and  physical  parameters  of  the  links  of  the  model;  2)  the  arbitrary  speci¬ 
fication  of  externally  applied  forces;  3)  the  use  of  relative  orientation 
angles  between  the  links  to  define  the  3vstera  configuration;  and  d)  the  use 
of  Lagrange' 3  form  of  d'Alembert's  principle  to  develop  the  governing 
equations  of  motion.  The  use  of  relative  orientation  angles  is  a  convenience 
in  the  specification  of  the  system's  configuration  and  in  the  introduction  of 


*  Numbers  in  brackets  refer  to  References  at  the  end  of  the  report. 


1 


flexible  and/or  torsional  springs  and  dampers  between  the  links. 

Lagrange's  form  of  d'Alembert's  principle  as  exposited  by  Kane  and 
others  [4-8],  has  been  shown,  for  large  systems,  to  possess  the 
advantages  of  both  Lagrange's  equations  and  Newton's  laws,  but  without 
the  corresponding  disadvantages.  That  is,  the  principle  provides  for 
the  automatic  elimination  of  the  "non-working"  internal  constraint 
forces  without  introducing  tedious  differentiation  or  other  similar 
calculations. 

The  objective  of  this  report  is  to  present  several  sets  of  data 
validating  this  cable  model.  This  data  consists  of:  1)  a  comparison  of 
results  obtained  from  this  model  with  analogous  results  obtained  from  a 
two-dimensional  multi-link  pendulum  model  with  governing  equations 
developed  by  Lagrange's  equations;  2)  a  comparison  of  data  from  the 
above  models  with  the  displacement  and  natural  frequencies  of  a  hanging 
cable  with  data  obtained  analytically  from  a  linear  partial  differential 
equation  model;  and  3)  a  comparison  of  model  data  for  a  submerged 
pendulum  with  experimental  data  recorded  at  the  Civil  Engineering  Laboratory 
at  Port  Hueneme,  California. 

The  balance  of  this  report  is  divided  into  four  parts  with  the  following 
part  summarizing  the  basic  equations  of  the  kinematics  and  dynamics  of  the 
cable  model.  This  is  followed  by  the  development  of  the  two-dimensional 
Lagrange  multi-link  pendulum  model.  The  comparisons  of  the  models  with 
each  other  and  with  analytical  and  experimental  data  are  presented  in  the 
next  part.  The  final  part  contains  some  concluding  remarks  on  the  signi¬ 
ficance  of  the  validation  and  on  the  application  of  the  model. 


THE  FINITE  SEGMENT  MODEL 


Configuration 

Consider  again  the  representation  of  the  model  as  shown  in  Figure  1. 
To  establish  a  geometrical  accounting  for  this  system,  select  one  link 
(say  the  first  link)  as  the  reference  link.  Next,  label  or  number  the 
other  links  in  ascending  progression  away  from  this  link  as  shown  in 
Figure  1.  The  configuration  and  kinematics  of  each  body  of  this  system 
may  then  be  developed  relative  to  the  reference  link  which,  in  turn,  has 
its  configuration  and  kinematics  defined  relative  to  an  inertial 
reference  frame  R. 


Consider  a  typical  pair  of  adjoining  links  such  as  L^  and  as  shown 
in  Figure  2.  Let  j<k,  that  is,  let  j»k-l.  Then  the  general  orientation 
of  L^  relative  to  L.. ,  may  be  defined  in  terms  of  the  relative  inclination 
of  the  dextral  orthogonal  unit  vector  sets,  n  and  n^  (i*l,2, 3)  fixed 
in  L_,  and  L^,  as  shown  in  Figure  2.  Specifically,  let  L^  and  L^  be  oriented 
so  that  n^  and  n^  are  respectively  parallel.  Then  L^  may  be  brought  into 
any  given  orientation  relative  to  L^  by  three  successive  dextral  rotations 
about  axes  parallel  to  o^,  o^.  and  a,  ^  through  the  angles  3,^,  ar.c  . 
n,.  and  n,  .  are  then  related  to  each  other  as: 

~jt  -tCl 


~ji 


SK 


im^km 


(1) 


3 


where  SK  is  a  3x3  orthogonal  transformation  matrix  ("shifter")  defined 
as  [9]  : 


SK 


im  -ij 


:m 


(2) 


(Regarding  notation,  repeated  subscripts,  such  as  m  in  the  right  side  of 
Equation  (1)  represent  a  sum  over  the  range  (1,2,3)  of  that  index.) 


SK  may  be  written  as  the  product  of  three  orthogonal  matrices  as: 


SK  =  AK  BK  GK 


where 
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(4) 


(5) 


(6) 


where  3  and  C  represent  the  sine  and  cosine. 
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These  expressions  allow  for  the  transformation  of  components  of 
vectors  referred  to  one  link  of  the  system  into  components  referred 
co  any  other  link  of  the  system,  and,  in  particular,  to  the  inertial 
reference  frame  R. 


Since  these  transformation  matrices  play  a  central  role  throughout  the 

analysis,  it  is  helpful  to  also  have  an  algorithm  for  this  derivative, 

expecially  the  derivative  of  SOK,  the  transformation  matrix  between  n,  ^ 

and  n  .  where  the  n  .  are  fixed  in  R.  Specifically,  SOK  is 
-ox  'Oi 


SOK 


ij 


n  . 


-01 


(7) 


Hence,  since  the  n  are  fixed  and  therefore  constant:  in  R,  the  following 
is  obtained: 

R 

d(SOK^)/dt  =*  nQi  •  dn^./dt  (8) 


where  the  superscript  R  indicates  that  the  derivative  is  computed  in  R. 
However,  since  n^  are  fixed  in  L^,  their  derivative  may  be  written  as 
x  where  is  the  angular  velocity  of  L,^  in  R.  Equation  (8)  may 
then  be  rewritten  as : 


dCSOK. . ) / dt 
ij 


*  j  n. 
'imn  kn-om 


-ki 


(9) 


or  as 


d ( SOK) / dt  *  WK  SOK 


(10) 


where  WK  is  a  matrix  defined  as 


VK 


im 


-e.  a. 
imn  kn 


(11) 


where  w,  are  the  n  components  of  w,  and  e.  is  the  standard  oermuta- 
kn  -on  -k  imn 

tion  symbol  [9,10].  WK  is  simply  the  matrix  whose  dual  vector  [10] 
is  Equation  (10)  thus  shows  that  the  derivative  may  be  computed  by 

a  simple  matrix  multiplication. 

Finally,  Figure  2.  contains  symbols  not  yet  defined  in  the  sequel. 

0_.  (j*l,...,N)  is  the  reference  point  of  ,  and  it  is  the  common  point 

of  L.  and  its  adjacent  lower-numbered  link.  G^.  (j*l,...,N)  represents 

the  mass  center  of  L j .  r^  (j=l,...,N)  is  the  position  vector  of 

relative  to  0.,  and  C,  is  the  position  vector  of  0,  relative  to  0,.  r. 

1  -k  k  j 

and  C,  are  thus  fixed  in  L.. 

-k  j 

Kinematics 


The  system  shown  in  Figure  1.  will  in  general  have  3N+3  degrees  of 
freedom.  These  may  be  defined  in  terms  of  generalized  coordinates  X 
(j*l , . . . , 3N+3) .  X^,  X^ ,  and  X^  represent  the  position  coordinates  of 

0.  in  R.  The  succeeding  triplets  of  coordinates  represent  dextral  rotations 
of  the  links  relative  to  the  respective  adjacent  lower  links. 

If  the  axial  rotation  of  the  links  is  neglected,  the  number  of 
degrees  of  freedom  can  be  reduced  to  2N+3.  This  also  avoids  singularities 
which  are  occasionally  encountered  with  large  rotations  as  discussed  in 
References  [11]  and  [12] .  When  the  axial  rotations  are  neglected, 


3 


Che  3,  are  zero  for  each  link  and  the  3K  become  idencitv  matrices.  In 
k 

che  following  kinematic  analysis,  the  axial  rotations  of  the  links  are 
neglected. 

The  angular  velocity  ui.  of  L,  in  R  is  readily  obtained  from  the 

""  !£  £ 

addition  formula  [6] : 


?2 


+  u. 


(12) 


where  u,  is  the 
-k 

transformation  properties  of  the  shifters,  may  be  writ 


angular  velocity  of  L,  relative  to 


By  using  the 
ten  as  : 


?k 


SOJ  (ad  + 
im  k  mi 


Y  AK  -)n 
k  m3  -oi 


(13) 


where  5  is  Kronecker's  delta  symbol  or  che  identicv  tensor  [9,10]. 
tnn 

Hence,  by  repeatedly  substituting  from  Equation  (13)  into  Equation  (12), 

ui,  takes  the  form: 

-k 


?k 


ai,  „  X,n 
k£m  i-ora 


(14) 


where  there  is  a  sum  from  1  to  2N+3  on  l  and  from  1  to  3  on  m.  From 
Equation  (13),  it  is  seen  that  the  non-zero  take  one  of  the  two 
forms : 


JkZm 


SOJ 


ml 


SOJ  AK  , 
mn  n3 


(15) 


depending  upon  whether  X,  is  a,  or  y,  . 


The  angular  acceleration  of  in  X  may  be  obtained  by  differ¬ 


entiating  Equation  (14).  Noting  that  the  n  are  constant,  this  becomes 


•  • 


3,  *  (w.  X,  +  u  X,)n 
~k  kirn  l  kirn  l  'om 


(16) 


where  from  Equation  (15)  the  non-zero  w,  ,  take  one  of  the  two  forms: 

Kim 


SOJ 


ml 


k£m 


(17) 


SOJ  AK  ,  +  SOJ  AK  . 
mn  n3  mn  n3 


where  SOJ  is  given  by  Equation  (10)  and  where  AK  is  obtained  by  differ¬ 
entiating  Equation  (4) . 


The  velocity  Vj  of  G_.  in  R  may  be  obtained  by  differentiating  , 


the  position  vector  of  G.  relative  to  a  fixed  ooint  0  in  R.  From 

1 


Figure  1.,  may  be  expressed  as: 


j-1 


h  •  (\  +  S0JkicJ1  * 


(18) 


where,  as  before,  there  is  a  sum  over  k  and  l  from  1  to  3.  Hence, 


V.  may  be  written  in  the  form: 
-J 


V.  =  V.  X.n  , 

-j  jik  L -OK 


(19) 


where  by  Equations  (13),  (10),  and  (11),  the  non-zero  V  are  given  by: 

J  *•* 


'  j Ik  “  ' ik 


(j  =  1,...,N;  l,k  =  1,2,3) 


(20) 
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¥ 


and 


Vjik  *  ';Jkp2.rjp  +  M-1  ^pi’Mp 


(21) 


where  WJ is  defined  as: 

£kp 


WJ,  ,  -  [3WJ.  /axJSOJ  (ft— 1, . . .  ,2N+3;  k,  P-1,2, 3)  (22) 

Kpl  kq  Z  qp 


Using  Equation  (11),  WJ^  ^  may  be  written  in  the  form: 


WJ,  -  -e,  SOJ 

kp  l  kqs  3  is  qp 


(23) 


The  acceleration  a^  of  in  R  may  be  obtained  by  differentiating 
Equation  (19),  leading  to: 


«  • 


?j  *  <W»  +  VW 


(24) 


where  by  Equations  (20)  and  (21) ,  the  non-zero  Vj  ik  are  given  by: 


*  '  i  •  M 

V.  -  WJ,  rJ  +2  WM, 

3  ik  kp  i  p  M=,1  kp^.  P 


(25) 


where  bv  Equation  (23)  WJ  is: 

kp  i- 


WJ  ,  *  -e,  (j).,  SOJ  +  oj.i?  SOJ  ) 

kp*  kqs  3  is  qp  3  *s  qp 


(26) 


Therefore,  the  kinematical  description  of  the  system  is  defined  by 
Equations  (Id) ,  (16),  (19),  and  (2d),  and  specifically  by  the  four  block 
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matrices  V^,  an<*  VjZk-  From  E<luacions  d5)>  d7)>  (21), 

and  (25),  it  is  seen  that  each  of  these  matrices  may  be  computed  by 
vector  and  matrix  multiplications  which  are  easily  developed  into 
computer  algorithms.  These  matrices  play  a  central  role  in  the  develop¬ 
ment  of  the  equations  of  motion  of  the  model. 


Equations  of  Motion 


Consider  again  the  cable  model  of  Figure  1.  Let  the  externally 
applied  force  system  on  each  link  be  replaced  by  an  equivalent  force 
system  consisting  of  a  single  force  F^ ,  passing  through  together  with 
a  couple  with  torque  .  Then  Lagrange's  form  of  d'Alembert's  principle 
states  that  the  governing  dynamical  equations  of  motion  for  the  chain 
system  are : 


Fz  +  F*  -  °  (£  -  1, . . . , 2N+3) 


F^  are  called  "generalized  active  forces"  and  they  are  given  by: 


F,  »  V  f  +  a)  M., 
I  3  ik  jk  j  ik  jk 


where  there  is  a  sum  from  1  to  N  on  j  and  from  1  to  3  on  k  .  and  where 
and  M.,  are  the  n  ,  components  of  F.  and  M. .  ?*  are  called 
"generalized  inertia  forces"  and  they  are  given  by: 


V . , ,  F*  +  jj  . . ,  M* 

J  ‘k  :k  jtlc  ].< 
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where  Che  indices  follow  the  same  rules  as  in  Equation  (28)  and  where 
and  are  the  components  of  the  inertia  forces  F*  and  inertia 
torques  M*  given  by  the  expressions  [6] : 


(no  sum) 


(30) 


and 


M*  »  x  (1^  *  cj^)  (no  sum) 


(31) 


where  m.  is  the  mass  of  L.  and  I  is  the  inertia  dyadic  of  L,  relative 
3  13  j 


to  G.  (j-l,...,N), 


By  substitutint  Equations  (16)  and  (24)  into  Equations  (30)  and 
(31)  and  ultimately  into  Equation  (27)  ,  the  equations  of  motion  may  be 
written  in  the  form: 


a„  X 
ZP  P 


( l  -  1, . . . , 2N+3) 


(32) 


where  there  is  a  sum  from  1  to  2N+3  on  p 
given  by : 


and  where 


a.  and  f,  are 

Z.O  i 


V 


+  I 


3  jpk  3  ^  3kn  3?n  j  ?'1<c 


(33) 


and 


■(F,  +  m,  V, V  .  X  +  I.,  a) 
L  f  Uk  ak  a  ikn 


•  • 


JZk’’jqk“q  '  *jkn“j  lnJ  jqkXq 


+  e  ,  j  .  jj,  -j  . , ,  I  XX) 
mnk  jqn  jsr  j  £ <c  jmr  q  s 


04) 


where  there  is  a  sum  from  1  to  N  on  j,  q,  and  s  and  a  sum  from  1  to 
3  on  the  other  repeated  indices. 


Equations  (32)  form  a  set  of  2N+3  simultaneous  ordinary,  nonlinear 

differential  equations  determining  the  2N+3  generalized  coordinates  Xj 

of  the  cable  system.  Since  the  coefficients  a,  and  f, ,  of  these 

4p  l 

equations  are  algebraic  functions  of  the  physical  parameters,  and  the 


four  arrays  V  V.  and  the  equations  may  be  generated 

on  a  computer.  Furthermore,  once  they  are  developed,  they  may  also  be 
solved  numerically  by  a  computer  by  using  a  standard  numerical  integra¬ 
tion  routine. 


Computer  Code 

Numerical  algorithms  to  evaluate  the  above  parameters  and  expressions 
have  been  developed  and  compiled  into  a  user-oriented  computer  code.  As 
input,  the  code  requires:  the  number  of  links;  the  masses;  the  centroidal 
principal  inertia  matrices;  the  mass  center  positions;  the  connection  (or 
reference)  point  positions;  the  motion  profile  for  those  links  with 
specified  motion;  the  applied  forces  and  moments:  and  the  initial 
configuration. 

The  output  of  the  code  includes:  the  values  of  all  variables  and 
their  first  derivatives;  the  mass  center  positions,  velocities,  and 
accelerations;  and  the  connection  point  positions,  velocities,  and 


accelerations. 


LAGRANGE  MULTI-LINK  PENDULUM  MODEL 


To  obtain  an  analytical  verification  of  the  governing  aquations  (32), 
consider  the  two-dimensional  oscillations  of  a  multilink  hanging  pendulum 
with  a  concentrated  end  mass  as  shown  in  Figure  3.  This  system  has  N 
degrees  of  freedom  which  may  be  described  by  the  angles  as 

shown  in  Figure  3.  The  equations  of  motion  of  this  relatively  simple 
system  can  be  obtained  through  Equations  (32)  or  independently  by  using 
Lagrange's  equations. 


If  each  link  has  the  same  mass  m  and  length  l  the  kinetic  energy 
K  of  the  system  may  be  expressed  as: 


K  -  <m/2)tv:  +  V2r  +  .  .  .  +  vl  +  (12/12)(?;  +  ?2  +  .  .  . 

G1  ~G2  g>j  12 

+  32)I  +  (M/2)V2  (35) 


where  the  (i  «  1,...,N)  are  the  mass  centers  of  the  links,  Q  is 
the  end  point  with  mass  M.  The  velocity  of  a  typical  mass  center  G 
may  be  expressed  as: 


Vi} 


■5) 


where  the  n^  (i  *  1,...,N)  a-e  unit  vectors  normal  to  the  links  as  shown 
in  Figure  3.  3y  substituting  expressions  of  the  form  of  Equation  (3h) 


inco  Equation  (35)  and  by  carrying  out  the  scalar  products,  the  kinetic 
energy  may  be  written  in  the  fora: 


K 


o  -V  M 

(ml1.1 2)  I  I  m,  . 

i-1  j-1  J 


where  the  m 


ij 


are  given  by  the  expressions: 


a  -  (1/2)  [1  +  2(S-k)  +  2(M/m)3cos(ej  -  1 ±) 

if  i#j  and  k  is  the  larger  of  i  and  j 


(37) 


(38) 


and 


m  -  C(N  -  i)  +  (1/3)  +  (M/m)  1  if  i-J  (39) 

Lagrange's  equations  then  lead  to  governing  equations  of  motion  of  the 
fora  [6]  : 


dOK/30  )dt  -  3K/30  -  F  (i  -  1,...,M) 

X  X  3  , 


(40) 


where  F  is  the  generalized  active  force  associated  with  3..  The  onlv 
9i  1 

forces  making  a  non-zero  contribution  to  Fa  are  the  weight  forces. 


Hence,  F,  may  be  written  as: 


vi 


'i 


>f 

Fa  -  z  (3V  /  304  )  •  (-mgk)  *  ( 3V  / 39  )  •  (-Mgk)  (41) 

'i  j-i  Mi  J  '  3 


where  k  is  the  vertical  unit  vector  aa  shown  in  Figure  3.  and  g  is  the 
gravitational  constant. 


3y  substituting  Equation  (37)  and  (41)  into  Equations  (40),  and 
by  carrying  out  the  indicated  operations,  the  equations  of  motion  take 
the  form: 


1'1  'Vj  *  VJ  +  • 0 


i  -  1 . N 


(42) 


where  the  are  given  in  Equations  (38)  and  (39)  and  where  p^  and 
k^  are  given  by  the  expressions: 


-  (1/2)  [l  +  2 (S-lcj  +  2(M/m)]sinOi  -  9 . ) 


(43) 


if  i  i  j  and  k  is  the  larger  of  i  and  j 


and 


pij  *  0  if  i  -  J 


(44) 


and 


ktj  -  0  if  i  *  i 


(45) 


and,  finally 


k  •  (M  -  i  -*•  (1/2)  (M/m) J sin9^ 


(46) 


MODEL  VALIDATION  AND  COMPARISON 


It  is  relatively  easy  to  show  -  particularly  if  the  number  of  links 
is  small  -  that  the  equations  of  motion  given  by  Equations  (32)  and 
(42)  are  the  same.  (This,  of  course,  requires  the  conversion  of  the 
relative  orientation  angles  used  in  Equations  (32)  to  absolute  angles  as 
shown  in  Figure  3.)  Hence,  when  Equations  (32)  and  (42)  were  independently 
integrated  for  a  number  of  cases  as  described  below  the  results  were 
identical . 

The  Hanging  Chain 

The  small  oscillations  and  natural  frequencies  of  a  hanging  cable  or 
chain  have  been  examined  and  studied  by  several  writers,  including 
Salvador!  and  Schwartz  [13]  and  Woodward  [14] .  These  investigations 
involve  the  solution  of  the  governing  partial  differential  equation 
modelling  the  continuum  of  the  hanging  cable. 

For  an  initially  straight  cable  inclined  at  an  angle  with  the 
vertical,  and  released  from  rest,  the  horizontal  displacement  y  may  be 
expressed  approximately  as  [13] : 

y  3  v(x,t)  3  80qL[O.139Jq(2.4O»'x/L)cos  1.20vg7L  t 

-  0.0175  Jq (5 . 32*'x/L)cos  2.76/g/'L  t  (47) 

+  0.00568Jq(3.65v’x/L)cos  4.33*g/L  t] 
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where  x  is  the  vertical  coordinate  along  the  cable,  t  is  the  time,  L 
is  the  total  cable  length,  and  Jq  is  the  Bessel  function  of  order  zero. 

For  a  15  ft.  cable  or  chain  modelled  by  15,  1  ft.  links,  Equation 
(32)  was  integrated  for  various  0^  and  the  results  were  compared  with 
those  produced  by  Equation  (47).  The  comparisons  of  the  predicted 
shapes  as  the  chain  passes  the  vertical  are  shown  in  Figures  4a.  and  4b. 

Matural  Frequencies 

For  small  oscillations.  Woodward  Cl4j  has  solved  the  governing 
partial  differencial  equation  of  a  hanging  cable  with  a  concentrated 
end  mass.  He  has  calculated  and  tabulated  the  natural  frequencies  for 
various  end-mass  to  cable-mass  ratios.  By  linearizing  Equation  (42), 
the  analogous  finite  segment  eigenvalue  problem  can  be  solved  and  the 
results  compared  with  those  of  Woodward  [14],  Table  I.  shows  a 
comparison  of  results  for  the  lowest  frequency  for  chains  of  5,  10,  and 
100  links  respectively.  Table  II.  shows  the  same  comparison  for  the 
second  lowest  frequency.  Finally,  to  show  the  convergence  of  the  finite 
segment  model,  Table  III.  presents  a  comparison  for  the  fifth  frequency 
for  chains  of  5,  10,  20,  50,  and  100  links.  (The  numbers  listed  in 
these  tables  are  2jj»L/g  where  j  is  the  natural  frequency.) 
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Knd-Mass/Cable-Muss  Ratio 


a3L2  I. 


A  Compar Ison  -f  the  Finite  Segment  Model  and  a  Ccncinuu"' 
Model  for  Che  lowe.sc  Natural  Frequery./  lor  a  "angiug  Cable. 
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_ 
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0.5 

2.174 
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1.0 
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2.113 
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j 
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2.067 
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I  5.0 

i 
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_ _ _ ; 

1 

1  10. 0 

— 

2.016 

2.016 

2.016 

_ 

2.016  1 

1  20.0 
i 

2.008 

2.008 

2.011 

_ _ — - 

\  ~ 

j  50.0 

i 

2.003. 

2.003 

r 

2.003 

2.001  : 

■ 

Natural  Frequence:  2~,  *"L/g 
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End-Mass/Cable-Mass  Ratio 


TABLE  II. 


A  Comparison  of  Che  Finite  Segmenc  Model  and  a  Continuum 
Model  for  che  Second  Lowest  Natural  Frequency  for  a  Hanging 
Cable. 


\.  N 

Number  of  Links 

Continuum 

M/mV 
c  V 

5 

10 

100 

Results  [14] 

0.0 

5.714 

5.574 

5.521 

5.520 

0.1 

5.788 

5.672  ' 

5.632 

5.632 

0.2 

6.067 

5.963 

5.929 

5.928 

0.5 

6.972 

6.873 

6.840 

6.840 

1.0 

8.308 

8.200 

8.164 

8.163 

2.0 

10.485 

-  - 

10.354 

10.311 

10.311 

5.0 

15.247 

15.060 

14.997 

1 

14.998  | 

10.0 

20.892 

20.637 

20.564 

20.553 

20.0 

29.059 

28.706 

28.614 

28.588 

50.0 

i 

45.478 

44.937 

44.673 

44.739 

Natural  Frequency  Zu^'L/'  g 


on 


TABLE  III.  A  Comparison  of  the  Finite-Segment  Model  and  a  Continuum  Model 
for  the  Fifth  Natural  Frequency  for  a  Hanging  Cable. 


Number  of  Links 


50 

■  Continuum  1 

Results  [14] | 

100  j  1 

i  ! 

15.065 

j 

14.968  |  14.931 

t 

17.686 

17.636  17.619 

.9.738 


19.722 


30.719 

26.223 

i  24.931 

! 

24.564 

24.512  | 

24.493  ! 

1 

37.291 

32.567 

31.014 

30.578 

30.515  | 

30.493  i 

47.904 

42.300 

40.311 

39.756 

39.675  1 

i 

39.648  | 

.  ■  ■  1 

70.871  !  62.869  59.928  59.112  58.996  58.951 

i  I 


97.861  112.282  |  82.820  j  81.679  81.546  j  81.466 


136.722  I  121.384  i  115.709  j  114.146  |  113.955  I  113.311 

i  i 


214.691  190.413  181.655  ! 179.264  | 173.654  178.511 


Natural  Frequency  2vj^»L/g 
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Submerged  Catenary  Cable 


An  experimental  verification  of  the  finite  segment  model  can  be  ob¬ 
tained  by  comparison  with  data  recorded  at  the  U.S.  Mavy  Civil  Engineering 
Laboratory  at  Port  Hueneme,  California.  In  these  experiments  a  totally 
submerged  cable  supported  at  one  end,  with  a  spherical  body  at  the  other 
end  was  initially  held  in  a  catenary  shape  and  then  released  from  rest 
[15,16].  The  subsequent  cable  shape  and  motion  were  recorded.  The 
same  experiment  was  simulated  using  the  finite  segment  cable  model  and  the 
computer  model  SEADYN  of  Reference  [15].  Table  IV  and  Figure  5.  show  a 
comparison  of  the  experimental  and  computed  results  for  the  end-body  displace¬ 
ment.  (In  the  analysis  using  Equations  (32),  the  fluid  forces  were  modelled 
by  using  results  of  References  [3]  and  [17].) 


■'ll  I  . . .  ■  I  MI  Him———'  -  — f  ;  .1  . J  ZEKr 


TABLE  IV.  A  Comparison  of  che  Finite  Segment  Model,  SEADYN  [15,16] 
and  Experimental  Results  for  the  End  Displacement  of  a 
Submerged  Cable. 


Time 

Experimental 
Results 
(X,Y) (ft) 

SEADYN 
Results 
(X,Y) (ft) 

v  mite-Segment 
Model  Results 
(X,Y) (ft) 

0.0 

20.5,  1.25 

20.5,  1.25 

20.5,  1.25 

2.0 

20.3,  6.65 

20.71,  8.86 

20.52,  7.97 

4.0 

20.2,  15.65 

21.01,16.81 

21.11,15.34 

6.0 

20.2,  24.15 

21.31,22.41 

21.0,  22.21 

8.0 

20.5,  32.65 

21.07,27.47 

20.89,28.97 

10.0 

20.5,  40.55 

20.77,32.29 

20.84,35.78 

_  . 

12.0 

20.5,  46.95 

20.54,37.35 

... 

20.38,42.61 

14.0 

20.5,  51.15 

20.24,42.11 

21.0,  49.19 

16.0 

20.0,  52.85 

20.24,46.81 

21.05,53.25 

18.0 

19.5,  53.55 

19.64,51.02 

20.63,54.06 

20.0 

18.9,  53.95 

19.17,53.07 

19.89,54.49 

_ 

NOTE:  X  is  the  Horizontal  Coordinate  of  the  Spherical 
End-Body  and  Y  is  the  Depth  Coordinate. 
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CONCLUSIONS  AND  APPLICATIONS 


Figure  4a.  shows  excellent  agreement  between  the  finite  segment 
model  and  the  linear  continuum  model  for  small  initial  displacement  angle. 
Indeed,  the  continuum  model  and  the  approximate  solution  of  Equation  (47) 
cannot  be  considered  valid  for  large  displacements.  Figures  4a.  and  4b. 
are  thus  in  a  sense  a  measure  of  the  range  of  validity  of  the  linear 
continuum  model.  Also,  note  that  for  large  9q  the  finite  segment  model 
shows  the  period  to  be  larger  than  that  predicted  by  the  linear  model. 

This  is  consistent  with  nonlinear  pendulum  theory  [18] . 

Tables  I,  II  and  III  show  the  convergence  of  the  finite-segment  model 
to  the  linear  continuum  model  for  small  amplitude  oscillations.  For  the 
lowest  frequency  there  is  excellent  agreement  with  only  five  links  in  the 
finite  segment  model.  Moreover,  even  for  the  fifth  frequency  there  is 
relatively  good  agreement  with  as  few  as  ten  links. 

Finally,  the  comparison  of  the  finite-segment  model  results  with 
experimental  results  also  shows  good  agreement.  This  comparison  validates 
not  only  the  finite-segment  model,  but  also  the  modelling  of  the  fluid 
forces  as  recorded  in  References  [3]  and  [17]  . 

These  results  all  suggest  that  the  finite-segment  cable  model  can 
provide  a  very  effective  and  efficient  model  of  the  nonlinear  dynamic 
behavior  of  long  heavy  cables.  Indeed  the  most  appropriate  applications 
are  likely  to  be  with  long  submerged  towing  cables,  mooring  cables,  and 


hoisting  cables. 
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